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ABSTRACT 

We report the detection of high-energy 7-rays from the quiescent Sun with the Large Area Telescope 
(LAT) on board the Fermi Gamma-Ray Space Telescope (Fermi) during the first 18 months of the mission. 
These observations correspond to the recent period of low solar activity when the emission induced by cos- 
mic rays is brightest. For the first time, the high statistical significance of the observations allows clear 
separation of the two components: the point-like emission from the solar disk due to cosmic ray cascades 
in the solar atmosphere, and extended emission from the inverse Compton scattering of cosmic ray elec- 
trons on solar photons in the heliosphere. The observed integral flux (>100 MeV) from the solar disk is 
(4.6 ± 0.2 [statistical error] j^J'g [systematic err or]) x 10 ^ cm ^ s ^, which is ^1 times higher than pre- 



dicted by the "nominal" modei of |Seckel et al. (1991). In contrast, the observed integral flux (>100 MeV) 
of the extended emission from a region of 20 radius centered on the Sun, but excluding the disk itself, 
(6.8 ± 0.7[stat.]iQ4[syst.]) x 10^^ cm^^ s^^, along with the observed spectrum and the angular profile, 
are in good agreement with the theoretical predictions for the inverse Compton emission. 
Subject headings: astroparticle physics — Sun: atmosphere — Sun: hehosphere — Sun: X-rays, gamma rays 
— cosmic rays — gamma rays: general 
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1. INTRODUCTION 



The Sun is a well-known source of X-ra^s an d 7-rays dur- 



ing solar flares (Peterson & Winckler 1959; Chupp et al. 1973 
[Kanbach et al. , 1 993 ), which are high-energy phenomena asso- 
elated with the flare-accelerated particle interactions in the so- 
lar atmosphere. Quiescent solar 7-ray emission from hadronic 
cosmic -ray (CR) interactions with t he solar atmosphere and 
photosphere was first mentioned by |Dolan & Fazio| ( |1965| l. 
[Pet erson et al. (1966) estimated its flux based on measure- 



ments of terrestrial emission, and Huds on (1989 ) suggested 
it to be detectable (^10""^ cm"^ s~^ above 100 MeV) by 
the EGRET experiment on board the Compton Gamma-Ray 
Observatory (CGRO). The first, and so far the only, detailed 
theoretical study of 7-ray emission from interactions of CR 
protons in the solar atmosphere was published by Seck el et al.| 
( 1991 1. The integral flux above 100 MeV was predicted to be 
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F{> 100 MeV) - (0.22 - 0.65) x 10"^ cm^^ g-i for their 
"nominal" model. However, attempts to observe such emis- 
sion with EGRET (1991-1995) yielded only an upper limit 
of 2.0 X 10"^ cm-^ s-^ a bove 100 MeV at 95% confidence 
level (Thom pson et al.|l"99 7). 

The existence of an additional, spatially extended compo- 
nent of the solar emission due to the inverse Compton (IC) 
scattering of CR electrons off sola r photon s was not real- 
ized u ntil recently ( Moskalenko et al. 2006 ; Orlando & Strong 
|2007| l. While the IC emission is brightest in the region within 
a few degrees of the Sun (hereafter we refer the angle relative 
to the Sun as the "elongation angle"), even at larger elongation 
angles it can be comparable in intensity t o the isotropic (pre - 
sumably extragalactic) 7-ray background ( Abdo et al. 2010c I. 
The flux for both components of the CR-induced emission is 
expected to change over the solar cycle due to the change of 
the heliospheric flux of the Galactic CRs in antic orrelation 
with the variations of the solar activity. Observations of the 
IC emission provide information about CR electron spectra 
throughout the entire inner heliosphere simultaneously, thus 
allowing comprehensive studies of the solar modulation of 
Galactic CRs in this region. For a moderately high - mod- 
erately low level of the solar modulation, the integral IC flux 
for el ongation angles <6° was predicted ( [Moskalenko et al" 



2006|l to beF(> 100 MeV) ~ (2.0- 4.3) 'x lO"'' cm'^ s 



respectively, and thus detectable by the EGRE T and Fermi- 
LAT. A calculation by |Orlando & Strong[ ( [2008 1) gave a similar 
flux F(>100 MeV) = 2.18 x 10"' cm for elongation 

angles <10° for the solar maximu m conditions. 

Reanalysis of the EGRET data by [Orlando & Strong [( |2008) 
led to the detection of both predicted components, point-like 
hadronic emission from the solar disk and extended leptonic 
emission from IC scattering of CR electrons on solar photons. 
Their analysis combined all 11 observational periods when 
the Sun was in the field of view of EGRET between the be- 
ginning of the CGRO mission (April 1991) and the end of the 
fourth observing cycle (October 1995). The average solar ac- 
tivity for this period was moderate, decreasing from its peak in 
1990 to the minimum in 1995. The [Orlando & Strong[ ( |2008]l 
analysis yielded fluxes F(>100 MeV) = (1.8 ± 1.1) xTTT^ 



cm-2 for the disk and F(>100 MeV) = (3.8 ± 2.1) x 
10^^ cm~^ s^^ for the IC component for elongation angles 
<10°, consistent with their estimate of the IC flux for the so- 
lar maximum conditions. 

The launch of Fermi in 2008 has made observations of the 
quiet Sun with high statistical significance and on a daily ba- 
sis possible. During the first two years of the Fermi mission, 
the solar activity has been extremely low, resulting in a high 
heliospheric flux of Galactic CRs. Therefore, the CR-induced 
quiescent 7-ray emission from the Sun is expected to be near 
its maximum. Preliminary analysis of the data from the Fermi 
Large Area Telescope (Fermi-LAT) showed the existence of 
both p oint-like and extended components of solar 7- ray emis- 
sion ( |Brigida||2009| |Giglietto|[2009l |Orlando[ [2009l). In this 
paper, we report on the Fermi-LAi observations of the qui- 
escent Sun during the first 18 months of the science phase of 
the mission. 



2. DATA SELECTION AND COUNT MAPS 

Fermi was launched on June 11, 2008 into circular Earth 
orbit with an altitude of 565 km and inclination of 25.6°, and 
an orbital period of 96 minutes. The principal instrument on 
Fermi is the LAT ( [Atwood et al.|2009j , a pair-production tele- 
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scope with a large effective area (^8000 cm^^ at 1 GeV) and 
field of view (2.4 sr), sensitive to 7-rays between 20 MeV and 
>300 GeV. After the commissioning phase, devoted to fine- 
tuning of the instrument and calibrations, the Ferm/-LAT be- 
gan routine science operations on August 4, 2008. The Fermi- 
LAT normally operates in sky-survey mode where the whole 
sky is observed every 3 hours (or 2 orbits) with an almost- 
uniform exposure on daily time scales. 

The energy dependent systematic uncertainties of the ef- 
fective area of the instrument were evaluated by comparing 
the efficiencies of analysis cuts for data and simulations of 
observations of pulsars (Rando 2009). This study revealed a 
systematic uncertainty of 10% at 100 MeV, decreasing to 5% 
at 560 MeV, and increasing to 20% at 10 GeV and above. 
The photon angular resolution is also energy dependent. The 
68% containment angle averaged over the Fenni-LPX ac- 
ceptance (the width of the point-spread function - PSF) can 
be approximated by the following expression: (0g8(e)) = 

([0.8°e"°-*]2 + [0.07°]2)^^^ where e is the photon energy 
in GeV. More details on the instrument performance can be 



found in the Fenni-hPX calibration paper (Atwood et al. 
2009[ l. The analysis presented here uses post-launch P6V3 



mstrument response functions (IRFs). These take into ac- 
count pile-up and accidental coincidence effects in the detec- 
tor subsystems that were not considered in the definition of 
the prelaunch IRFs. 

We use the Fermi-hP^ data collected between August 4, 
2008, and February 4, 2010. Events >100 MeV arriving with 
elongation angles < 20° (region of interest - ROI) an d sat- 
isfying the Diffuse class selection (Atwood et al. 2009) are 
used. To reduce the contamination by the 7-ray emission 
coming from CR interactions in the Earth's upper atmosphere 
our selection is refined by selecting events with zenith angles 
<105°. To reduce systematic uncertainties due to the bright 
diffuse 7-ray emission from the Galactic plane and a possible 
spillover due to the broad PSF at low energies, we have also 
excluded the data taken when the Sun was within 30° of the 
plane (I60I > 30°). We further excluded the periods when 
the Sun was within 20° of the Moon or any other bright ce- 
lestial source with the integral flux Fifgl > 2 x 10^^ cm^^ 
s^^ above 100 MeV as selected fr om the IFGL Fermi-LAiT 
source catalog ( |Abdo et al.|2010^ . These various selections 
produce a very clean event subsample but at the expense of re- 
moving about 93% of the initial ROI data set as summarized 
in Tabled 

Because the Sun is moving across the sky, the analysis of 
its emission requires special treatment. Therefore, a dedicated 
set of tools was developed, not a part of the standard Fermi- 
LAT Science Tools package, to deal with moving sources such 
as the Sun and the Moon. Using these specialized tools, the 
data are selected in a moving frame centered on the instanta- 
neous solar position, which is computed using an interface to 
the JPL ephemeris libraries^^. 



3. ANALYSIS METHOD 

For the analysis of the Sun-centered maps we used the 
Fermi-LAT Science Tools^** version 9rl6p0. The gtlike 
tool provides maximum likelih ood p a ramet er values (using 
the method described in CashlfT979l IStrong 1985; Mattox 



|et al.|1996l l, which derives error estimates (and a full covari- 
ance matrix) from Minuit, a minimization tool supported by 
CERN, using the quadratic approximation around the best fit. 

3.1. Background determination 

The correct evaluation of the background in the region 
around the Sun is of considerable importance for the analysis 
of the weak extended IC emission. The l atter is expected to 
decrease as ^\/9 with elongat ion angle (Moskalenko et al. 



http://iau-comm4.jpl.nasa.gov/access2ephs.html 

Available from Fermi Science Support Center (FSSC), 
[http://fermi.gsfc.nasa.gov/ssc 



|2006[ [Orlando & Strong '2007) and becomes indistinguish 
able from the background for 6 > 20°. The background is 
mainly due to the diffuse Galactic and isotropic (presumably 
extragalactic) 7-ray emission averaged along the ecliptic and 
to weak point sources. The evaluation of the background was 
done by two methods, one based on the analysis of flight data 
and the other on simulations. 

The first method utilizes the data and is called a "fake-Sun" 
analysis, where an imaginary source trails the Sun along the 
ecliptic. For this method, application of exactly the same 
sets of cuts as are used for the Sun yields an estimate of 
the background. Since the extended solar IC emission is in- 
significant for elongation angles 9 >20° we used 40° as the 
minimum trailing distance. To reduce statistical errors in the 
background determination, the background is averaged over 
4 fake-Sun sources displaced from each other and from the 
Sun itself by 40° intervals along the ecliptic. Because all 
fake-Sun sources are sampling the same area on the sky, the 
backgrounds determined using individual fake-Sun sources 
are consistent within a fraction of a percent. Note that the 18- 
month analysis period is long enough to average out any ef- 
fects connected with incomplete sampling of the background. 
This would be an issue if the analysis period was shorter than 
12 months. 

Figure [T] shows the Gaussian-smoothed count maps >100 
MeV centered on the solar position and the hypothetical trail- 
ing source (average of the 4 fake-Suns). The solar emission 
is clearly seen on the left panel, while the right panel shows 
the background, which is essentially uniform. The integral 
intensity distribution for the two samples, centered at the so- 
lar position and centered on the averaged fake-Sun source, is 
shown in Figure [2] The number of events per solid angle is 
shown vs. the angular distance from the Sun (the elongation 
angle) and the fake-Sun positions for a bin size 0.25°. While 
for the solar centered data set the integral intensity increases 
considerably for small elongation angles, the averaged fake- 
Sun profile is flat. The two distributions overlap at distances 
larger than 20° where the signal significance is diminished. 
The gradual increase in the integral intensity for >25° is 
due to the bright Galactic plane broadened by the PSF, see the 
event selection cuts summarized in Section |2] and Table [T] 

The second method of evaluating the background uses an 
all-sky simulation which takes into account a model of the 
diffuse emission (including the Galactic and isotropic compo- 
nents, gll jem_v02.fits and isotropic_iem_v02.txt, correspond- 
ingly; see footnote [58]l and the sources from IFGL Fermi- 
LAT catalog ( A bdo et al.|2010a 1). To the simulated sample we 
apply the same set of cuts as applied to the real data and se- 
lect a subsample centered on the position of the real Sun. The 
simulated background is then compared with the background 
derived from a fit to the fake-Sun in the first method. Figure|3] 
shows the spectra of the background derived by the two meth- 
ods. The agreement between the two methods (and the spec- 
trum of the diffuse emission at medium and high latitudes, 
I Abdo et aL]2010c. not shown) is very good, showing that the 
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background estimation is well understood and that there is no 
unaccounted or missing emission component in the analysis. 

Finally, we check the spatial uniformity of the background 
determined by the fake-Sun method. The ROI restricted by 
9 < 20° was divided into nested rings. We use 4 annular rings 
with radii 9 ^ 10°, 14°, 17.3°, and 20°, which were cho- 
sen to subtend approximately the same solid angle for each 
ring, and hence should contain approximately equal numbers 
of background photons if their distribution is spatially flat. 
The ring-by-ring background intensity variations were found 
to be less than 1%. Note that the background emission is con- 
siderab ly m ore intense than the expected IC component (see 
Section [J!2] i, and even small background variations across the 
ROI may affect the analysis results. To minimize these sys- 
tematic errors, we therefore using the ring method for the 
background evaluation. 

The evaluated spectrum of the background for 9 < 20° was 
fitted using the maximum likelihood method and the results 
were used to derive the simulated average photon count per 
pixel using the gtmodel routine from the Ferm/-LAT Science 
Tools. For each pixel of the predicted average photon count 
distribution, we generate a set of 100 random events assum- 
ing Poisson statistics and compare them with the observed 
number of photons in the corresponding pixel of the map cen- 
tered on the fake-Sun, as shown in Figure |4] The resulting 
distribution of the difference of counts between the observed 
and simulated photons per pixel is consistent with a normal 
distribution of mean 0.155 ± 0.002 and standard deviation 
2.337 ±0.002. 

3.2. Separation of the solar emission components 

The separation of the disk and extended components of the 
solar emission was done in a model-independent way. The 
Sun-centered maps of Fermi-hP^ 7-ray counts were ana- 
lyzed using a maximum likelihood technique where the back- 
ground parameters in each nested ring were fixed using the 
fake-Sun method as described above. The flux and spectral in- 
dex of the extended emission were determined independently 
in each energy range in each nested ring, while the disk com- 
ponent (modeled as the Fermi-LP^ PSF because the ^0.5° 
solar disk is not resolvable) was allowed to have a free spectral 
index and flux normalization. The choice of the annular radii 
(5°, 11°, and 20°) has been optimized to have a likelihood test 
statistic^*^ (TS) value > 25 for the fitted IC component in each 
ring. 

Figure [5] shows the angular distribution of photons 
>500 Me V from the Sun on a 0.2° grid, the background de- 
termined by the fake-Sun method, and the fitted disk and ex- 
tended components. The spectrum of a disk source is modeled 
as a power law with the total flux and spectral index obtained 
from the fit of the disk component. The extended emission is 
modeled as the sum of the individual fits in each energy bin 
in each nested ring. The observed angular distribution can be 
well fitted only by adding an extended IC component. 

4. RESULTS 

The analysis of the ferm/-LAT observations for the first 18 
months of the mission gives a significant detection and sep- 
aration of the two components of solar 7-ray emission. The 
large photon statistics allow us to derive the spectral shape 
of each component by fitting them in narrow energy bands. 



so that the exposures and convolutions with the PSF do not 
depend strongly on a the assumed spectral shape. 

Table [2] gives the energy spectrum for the IC component 
from the model-independent analysis or elongation angles 
<5°, which corresponds to the radius of the innermost ring, 
and <20°. The fitted integral fluxes for the IC component 
in each ring are given in TableJS] together with model cal- 
culations (described in Section pF The energy bin size was 
selected to provide good convergence of the likelihood fit in 
each energy interval. The observed spectrum for the disk 
component is given in Table [4] 

Evaluation of the systematic errors for each component of 
the solar emission has to take into account the uncertainties 
in the effective area as a function of energy as indicated in 
Section|2]and the statistical uncertainties in the determination 
of the background (Section [3T| i. The uncertainties associated 
with the effective area were propagated using the modified 
effective areas brack eting the nominal ones (P6V3 Diffuse) 
(Abdo et al. 2010b| l. The uncertainties in the background 
were taken into account by assigning ±.la deviation to the 
background flux and repeating the fit for all components. 

5. CALCULATION OF THE IC EMISSION 

The calculation of the IC emission from the solar r adia tion 
field was first de scribed by Moskalenko et al. ( 2006| l and Or- 
lando & StrongI (2008) using the formula for the differential 
interaction rate for an anisotropic distribution of target pho- 
tons ( [Moskalenko & Strong|2000[ ). The CR electron distribu- 
tion was assumed to be isotropic (hereafter we use the term 
electrons to refer to both electrons and positrons). The up- 
scattering of optical solar photons to the Ferm;-LAT energy 
range >100 MeV involves CR electrons above ^2 GeV. At 
energies below ^20 GeV, CR electrons in the heliosphere are 
subject to significant adiabatic energy losses and drifts in the 
magnetic field: th e combined effe ct is called the heliospheric 
modulation (e.g., |Potgieter]|1998| l. A calculation of the spa- 
tially extended ernission due to the IC scattering of CR elec- 
trons on solar photons requires the integration of 7-ray yields 
along the line-of-sight folded with the CR electron spectru m 
at different heliospheric distances ( (Moskalenko et al.|2006|l. 

Due to the radial distribution of the solar photons the IC 
emission is brightest at small elongation angles, and is dis- 
tributed over the whole sky at low levels. The differential flux 
of upscattered photons depends on how the CR electron spec- 
trum changes with heliocentric distance, which is the only un- 
known in the calculations. Therefore, observations of the IC 
component of the solar emission provide a new tool to probe 
the CR electron spectrum in the inner heliosphere (r < 10 
AU) down to close proximity to the Sun. 

5.1. CR modulation in the heliosphere 
Studies of the CR modulation are based on the solution of 



The likelihood test statistic TS is defined in Mattox et al. ' 1996 '. 



the,Parker,(1965) transport equ ation (e.g., see reviews by Pot- 
gieter|1998y Heber et al.|2006| . Particle transport to the inner 
heliosphere is mainly determined by spatial diffusion, con- 
vection by the solar wind, and adiabatic cooling. Besides, 
rotation of the Sun causes the interplanetary magnetic field 
in the solar equatorial plane to be distinctly different from 
the fie ld above and below the poles, the so-called Parker's 
spiral (Parker 1958 1. This, in turn, causes differences in 
the spectra of CR particles in these regions. Realistic time- 
dependent three-dimensional hydrodynamic mod els incorpo- 
rating these effects have been developed (e.g., ,Florinski et aL| 
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2003] ILangner et al. 12006} [Potgieter & Langner|2004| i; how- 
ever, the effect of heliospheric modulation is still far from be- 
ing fully understood. One of the major difficulties in devel- 
oping models of heliospheric modulation is that the data gath- 
ered by spacecraft taken at different heliospheric distances are 
often at different energies and related to different modulation 
levels. The region <1 AU is the least studied. Another prob- 
lem is that the input information for the modulation models, 
such as the local interstellar spectra of CR species, is missing. 

Although much effort has gone into development of realis- 
tic modulation models, the method most often used is the so- 
called "force-field" approximation (Gleeson & Axford 1968|, 
which employs a single parameter - the "modulation poten- 
tial" $ - that varies over the solar cycle to characterize the 
strength of the modulation effect on the CR spectra: 

(E^ - m^c^) 

Ur.E)^M^^E + <,[r]) ^ , (D 

where Je(oo, E + ^[r])\s the local interstellar electron spec- 
trum, E is the total electron energy, mgC^ is the electron rest 
mass, $(r) is the modulation potential, and r is the distance 
from the Sun. Even though this approximation is very crude 
and implies spherical symmetry for the heliosphere, it can be 
used as a convenient parameterization of the CR spectrum at 
different levels of solar activity. 

Exp ressions for "I>(r) were derived by Moskalenko et al. 
(|2006 ) based o n the radial dependence of the CR mean free 
path given by Fujii & McDonald] ( [2005] l for Cycles 20/22 and 
and for Cycle 21 assummg separability of the heliospheric dif- 
fusion coefficient into radial and energy-dep endent functions. 
For Cycles 20/22, [Moskalenko et al^ ( |2506) l obtain: 

i:88\0.24 + 8(r-o i-r„-° i), r < ro, 

where $o is the modulation potential at 1 AU, rg — 10 AU, 
and rf, — 100 AU is the heliospheric boundary. For Cycle 21 
they give: 



*2W = %{r' 



(3) 



These formulae were derived for r > 1 AU. Closer to the 
Sun at 0.3 AU ^ r < 1 AU there are only very few mea- 
surements of Galactic CR protons and heliu m, and these are 
at low energies <60 MeV nucleon"^ ( e.g., Christon et al. 
[1975) IKunow et all[T977l [Miiller-Memn et al.,,1977j . Smce 
CR transport is strongly energy dependent, these low-energy 
protons and helium nuclei measurements are irrelevant to the 
current study and CR transport is very uncertain. Therefore, 
as a first approximation we use eqs. ([3]) for the entire he- 
liosphere (Models 1, 2), from the solar surface to the helio- 
spheric boundary Rq < r < ri,, except for Model 3 described 
below. 

5.2. The CR electron spectrum 

Assuming that the CR propagation in the heliosphere is 
spherically symmetric in the first approximation, the CR elec- 
fron specfrum at different heliospheric distances is the only 
unknown in evaluating the solar IC flux. Comparison of the 
model calculations with observations provides a method to 
probe the CR elecfron spectra at different heliospheric dis- 
tances r < 10 AU. Because of the lack of CR measurements 
outside of the heliosphere, we have to rely on a comprehen- 
sive Galactic CR propagation model tuned to othe r CR data 
and diffuse 7-ra y emissi on, such as G ALPROP (Str ong & 
[Moskalenko! 1998, Moskalenko & Strong. 1 998 ,, Ptuskin et al. 



[20061 lAbdo eTaL]|2009a| [20T0cl [Vladimirov et al.|2010l l, or 
measurements of the flux at 1 AU. 

In this paper we use the CR electron spectru m recently 
measured b y the Ferm/-LAT ( |Abdo et a l. 2009b ; Ackermann] 
et al.|2010| l between 7 GeV and --1 TeV. At energies <7 GeV 



Je{0O,E) 



(4) 



we u se tEe Alpha Magne tic Specfro meter (AMS-01) lepton 
data ( [Alcaraz et al.|2000[ ) collected in its flight in June 1998, 
which are the most accurate to date, but were made in the 
previous solar cycle. Preliminary results by the Payload for 
Antimatter Matter Exploration and Light-nuclei Astrophysics 
(PAMELA, Adriani et al. 2011), in orbit since June 2006, in- 
dicate that the electron spectrum >7 GeV is consistent with 
the Ferm/-LAT data, but is more intense than the AMS-01 
data at lower energies. 

Unusually high intensities of Galactic CRs during the pe- 
riod of anomalously low solar activity were also reported by 
the Advanced Composition Explorer (ACE), which monitors 
the flux of CR species at 1 AU in a few 100 MeV/nucleon en- 
ergy ran ge. The measured in tensities of major species from 
C to Fe ( Mewaldt et al.|20 10') were each 20%-26% greater in 
late 2009 than in the 1997-1998 minimum and previous solar 
minima of the space age (1957-1997). While the value of $0 
is somewhat arbitrary and depends on the assumed local in- 
terstellar spectrum (but does not depend on the particle type), 
the unusually low solar activity during the observational pe- 
riod suggests that the modulation of Galactic CRs should be 
considerably weaker than at any other time. 

Therefore, the Fen«/-LAT CR electron spectrum was fitted 
using a parameterization for the local interstellar spectrum 
and assuming a relatively low value of the modulation poten- 
tial $0 = $(1 AU) = 400 MV (eq. (l)): 

a{b + c)-<^{E/b)-^,E <b 
a{E + c)-'^,E> b, 

where Je is in units m^^ s^^ sr^^ GeV~^, a ~ 160.24 m~^ 
s-i sr-i GeV2 "3, b = 7 GeV, c = -1.20 GeV, d = 3.03, 
and E is the total electron energy in GeV. Different values 
for <i>o yield somewhat different values for the other param- 
eters while providing the same fit to the Fermi-LAI electron 
spectrum. We found that the value of $0 given above pro- 
vides reasonable agreement with the observed IC flux (see 
Section [4]l; more detailed analysis of the electron specfrum 
in the inner heliosphere based on ferm/-LAT observations of 
the IC emission will be given in a forthcoming paper 

Since most of the IC emission is produced by elecfrons in 
the inner heliosphere r < 10 AU, eqs. Q, Q is a good 
approximation of the CR electron spectrum even though it 
may produce an unphysical spectrum beyond the heliospheric 
boundary at r?, = 100 AU. (In fact, for elongation angles 
9 < 20° most of the e mission is produced by e lectrons at 
r < 2 AU, see Figure 2 inlMoskalenko et al. 2006 ) However, 
the electron spectrum obtained in such a way is also close to 
the local interstellar spectrum that is calculated by GALPROP 
and used for modeling of the Galactic diffuse 7-ray emission 
dAbdo et al. 2009a, 2010c]l. 

Using this parameterization for the CR electron spectrum, 
we construct three models to compare with the data. To obtain 
the electron spectrum at an arbitrary heliospheric distance, r, 
Models 1 and 2 use parameterizations ^i.2{r) (with $0 = 
400 MV) given by eqs. ([2]), ([3]l, respectively. Model 3 uses 
parameterization eq. ([2]i as specified above, but assuming no 
additional modulationfor r < 1 AU, i.e., $(r < 1 AU) = 
$0- In this model the electron specfrum at <1 AU is the same 
as measured by the Fermi-LAI, which provides us with an 
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upper limit for CR electron flux closer to the Sun. 

Figurel6]shows the fit to the Ferm/-LAT electron spectrum. 
The AMSTOl lepton data were collected during the period 
of rising solar activity''", where the low modulation poten- 
tial adapted for the current solar minimum is not appropri- 
ate. These data are not fitted and shown for only illustration. 
To illustrate the effect of different parameterizations ^i,2{r), 
the Figure also shows the spectra for Models 1, 2 at r = 0.3 
AU. All three models of the electron spectrum are normalized 
to the same Ferm/-LAT data at 1 AU. Therefore the largest 
difference between the models is at low energies where the 
modulation effect is strongest. 



6. EMISSION FROM THE SOLAR DISK 

A detailed theoretical s tudy of 7-ray emissi on from the so- 
lar disk was published by |Seckel et al ( 1991| l. The predicted 
flux is very sensitive to the assumptions about CR transport 
in the inner heliosphere, the magnetic fields in the solar atmo- 
sphere, and the CR cascade development in the solar atmo- 
sphere. 



The jSeckel et al. ( 1991 ) calculation divides space into three 
regions: the interplanetary space, the corona, and the regions 
below the corona. The interpla netary magnet ic field (>2Rq) 
is taken to be in the form of a Parker ( 1958) ) spiral with the 
total field ^50 /iG near the Earth. Inferior to Earth's orbit the 
field is nearly radial at angle ^45°. CR propagation from the 
Earth's orbit to the bottom of the corona is treated as spheri- 
cally symmetrical diffusion with no absorption (i.e., neglect- 
ing CR interactions). The adopted diffusion coefficient has a 
linear dependence on energy, while the spatial radial depen- 
dence is a power-law with index 2. The C R spectrum at 1 
AU was adopted from Webber & Potgieter|(1989), which is 
a parameterization of the measu rements made^Funng balloon 
flights in 1976 and 1979 ( [Webber et al.|1987l l. 

Interior to the corona, the CR propagation and absorption 
(interactions) are treated simultaneously using a Monte Carlo 
code. The magnetic field configuration is chosen correspond- 
ing to a quiet Sun (Priest 1982), where the average field 
strengths are of the order of a few Gauss. However, the fields 
are non-uniform with the flux bundles (~10^ G, a few hun- 
dred kilometers across) located at the corners of convective 
cells and separated by thousands of kilometers. The chromo- 
sphere is assumed to be isothermal in hydrostatic equilibrium 
which gives an exponential density profile with scale height 
~1 15 km. The density profile b elow the photosphere is taken 
from Baker & Temesvary ( 1966j). The characteristic column 



density is ~40 g cm ^ at 500 km depth and 100 g cm at 
900 km. 

The cascades initiated by high-energy particles (>3 TeV) 
do not contribute much to the observed 7-ray flux and are 
neglected. The solar magnetic fields do not significantly af- 
fect their directionality until particle energies drop below ~ 10 
GeV - by that time most cascades are deep enough so that 
only a few low-energy photons will escape. On the other 
hand, the low-energy primaries produce cascades for which 
a significant amount of energy is reflected back to the solar 
surface (so-called mirrored showers). The typical low-energy 
cascade has less than a few interaction lengths of material to 
pass through before photons can escape. Therefore, such cas- 
cades act as if they evolve in moderately thin targets. 

See the charts of the computed tilt angle of the heliospheric current sheet 
at http;//wso. stanford.edu/ 



The actual yields are calculated by propagation of one- 
dimensional cascades through a slab. The photon yield in- 
cludes only the photons that make it through the slab. Sec- 
ondary electrons, positrons, and baryons exiting the slab are 
ignored even though they are likely to reenter the Sun. This 
may u nderestimate the ac tual 7-ray flux. 

The jSeckel et al.| ( fl99I] l calculations were made in two sce- 
narios: (i) neglecting the effects of interplanetary magnetic 
field on particle propagation and assuming the solar surface is 
fully absorbing (so-called "naive" model), and (ii) the "nom- 
inal" model, which includes all the assumptions about CR 
diffusion in the inner heliosphere and corona. The integral 
flux above 100 MeV was predicted to be F{> 100 MeV) ~ 
(0.22 - 0.65) X 10"'^ cm-2 g-i for the "nominal" model, 
where the range corresponds to the different assumptions 
about the CR cascade development: slant vs. more realistic 
mirrored showers (i.e., reflected back to the solar surface). 

7. DISCUSSION AND CONCLUSIONS 

Figures |7] and [8] show the results of the analysis of the ex- 
tended emission component together with the model predic- 
tions (Section |5]l. The plotted values are obtained from the 
IC flux in each ring as shown in Table [3] divided by the cor- 
responding solid angle. The model calculations are shown 
unbinned (curves) and binned with the same bin size as used 
for the data. Although the highest energy point 3 — 10 GeV 
shows some excess relative to the model predictions, this is 
difficult to explain from the model viewpoint since the effect 
of the solar modulation is decreasing at high energies thus 
making the model more accurate. Future analysis with larger 
statistics should clarify if this discrepancy is real. The agree- 
ment of the observed spectrum and the angular profile of the 
IC emission with the model predictions (as described in Sec- 
tion |5]) below a few GeV is very good. The innermost ring 
used for the analysis of the IC emission subtends an angular 
radius of 5° corresponding to a distance ^0.1 AU from the 
Sun, i.e. 4 times closer to the Sun than Mercury. At such a 
close proximity to the Sun, and actually anywhere <1 AU, 
the spectrum of CR electrons has never been measured. 

It does not seem possible to discriminate between the mod- 
els at the current stage. The spectral shape <1 GeV in Fig- 
ure |7] and the intensity in the innermost ring in Figure [8] is 
better reproduced by Models 1, 2, while the intensity in the 
middle ring 5° — 11° (Figure[8]l is better reproduced by Model 
3. Even though the current data do not allow us to discrim- 
inate between different models of the CR electron spectrum 
at close proximity to the Sun, the described analysis demon- 
strates how the method would work once the data become 
more accurate. In particular, it is possible to increase the 
statistics by four-fold by masking out the background sources 
or modeling them, instead of requiring the angular separa- 
tion between bright sources and the Sun to be >20° (Table[T]). 
More details will be given in a forthcoming paper The in- 
crease of the solar activity may also present a better oppor- 
tunity to distinguish between the models since the difference 
between the model spectra of CR electrons will increase with 
solar modulation. 

The intensity of the IC component is comparable to the in- 
tensity of the isotropic 7-ray background even for relatively 
large elongation angles (Table [2]). Integrated for subtended 
angles <5 °, the latter yields ~ 2! 5 x 10^^ cm^'^ s^^ above 
100 MeV ( |Abdo et al.|20I0c| ) vs. -1.4 x lO"'^ cm^^ g-i for 
the IC component. For subtended angles <20°, the integral 
flux of the isotropic 7-ray background is ~3.9 x 10~^ cm~^ 
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above 100 MeV vs. -6.8 x IQ-^ cm" 



for the IC 



component. Therefore, it is important to take into account the 
broad nonuniform IC component of the solar emission when 
dealing with weak sources near the ecliptic. The relative im- 
portance of the IC component will increase with time since the 
upper limit on the truly diffuse extragalactic emission could 
be lowered in future as more 7-ray sources are discovered and 
removed from the analysis. 

Figure [9] shows the spectrum for the disk component mea- 
sured by the Fenni-LAT (Table p]) and two mod el predictions 
("naive" and "nominal") by Secke l et al.| ( 1991]l as described 
in Section [6] In each set of curves, the lower bound (dotted 
line) is the CR induced 7-ray flux for the slant depth model, 
and the upper bound (solid line) is the 7-ray flux assuming 
showers are mirrored (as charged particles would be). The 
observed spectrum can be well fitted by a single power-law 
with a spectral index of 2.11 ± 0.73. The integral flux of the 
disk component is about a factor of 7 higher than predicted 
by the "nominal" model. An obvious reason for the discrep- 
ancy could be the conditions of the unusually deep solar min- 
imum during the reported observations. However, this alone 
can not account for such a large factor, see a comparison with 
the EGRET data below. Another possibility for an estimated 
"nominal" flux to be so low compared to the FermZ-LAT ob- 
servations is that the secondary particles produced by CR cas- 
cades exiting the atmospheric slab are ignored in the calcu- 
lation while they are likely to reenter the Sun. On the other 
hand, the proton spectrum by Webber et al. ( 1987 ) used in the 
calculation is about a factor of 1 .5 higher above —6 Ge V than 
that measured by the B ESS experiment in 1998 (see Figure 
4 in |Sanuki et al.|[2000l ). Meanwhile, calculation of the disk 
emission relies on assumptions about CR transport in the in- 
ner heliosphere and in the immediate vicinity of the Sun thus 
allowing for a broad range of models (cf. "naive" vs. "nom- 
inal" models). The accurate measurements of the disk spec- 
trum by the FermZ-LAT thus warrant a new evaluation of the 
CR cascade development in the solar atmosphere. 

The spectral shape of the observed disk spectrum is close to 
the predictions except below —230 MeV where the predicted 
spectral flattening is not confirmed by the observations. This 
may be due to the broad PSF making it difficult to distinguish 
between the components of the emission or a larger systematic 
error below —200 MeV associated with the IRFs. 

The results of Ferm/-LAT observations can be also c om- 
pared with those from the analysis of the EGRET data ( |Or-| 
|lando & Strong! [2OO 8). The latter gives an integral flux 
(>100 MeV) for the disk component of (1.8 ± 1.1) x lO"'^ 
cm"2 s"^ for the IC flux Orlando & Strong (2008) obtain 
(3.9 ± 2.2) X 10"'^ cm-2 ^-i ^j. elongation angles <10°. 
The EGRET-derived integral IC flux is very close to the flux 
observed by the Ferm/-LAT for the same range of elonga- 
tion angles (Table|3]l. Meanwhile, the integral flux of the disk 
component observed by the Fermi-LNI (Table |4|i is a factor 
of —2.5 higher than obtained from EGRET data. Such an 
increase in the 7-ray flux from the solar disk may imply a 
significant variation of the disk emission over the solar cycle. 
This is not surprising since the disk flux depends on the am- 
bient CR proton spectrum in the immediate proximity of the 
solar atmosphere, in contrast to the solar IC emission, which 
is produced by electrons in a considerable part of the helio- 
sphere and integrated along the line of sight. However, such 
variations should be confirmed by future observations. 

The analysis reported here refers to the period of the so- 



lar minimum. However, the solar activity is just starting to 
pick up for the current cycle. We expect that the effects of 
the increased solar activity on the fluxes of the two compo- 
nents of the emission will be different and it is important to 
follow their evolution during the whole solar cycle. Continu- 
ous monitoring of the solar emission by the Ferm/-LAT over 
the whole solar cycle will enable us to study CR transport in 
the inner heliosphere, to improve models of the solar mod- 
ulation, and of the development of CR cascades in the solar 
atmosphere. 
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Table 1 

Summary of the event selection cuts 



Cumulative event selections 




Photons, 

% 


Livetime, 

% 


eRoi < 20° 




100 


100 


\bQ\ > 30° ^ 




29.2 


60.9 


^Moon > 20° 




26.2 


54.8 


&{FiYi,L > 2 X 10 " cm - s 


;^ 20-' 


6.5 


17.6 



Table 2 

Differential spectrum of the IC component 



Energy range, 


Fluxitstatitsyst, cm s 


1 MeV-i 


(E^Fluxjitstatitsyst, 10"^ cm-^ s"! MeV 


IQS MeV 


61 < 5° 


e < 20° 


6» < 5° 6> < 20° 



0.1 - 


0.3 


(4.4 ± 0.9 ± 0.1) X 10-1" (2.1 ± 0.5 ± 0.3) x 10 


-9 


1.3 i 0.2 it 0.03 


6.2 it 1.4 it 0.1 


0.3 - 


1.0 


(4.9 ± 1.0 ± 0.1) X 10-" (2.2 ± 0.5 ± 0.3) x 10" 


-10 


1.4 it 0.3 it 0.03 


6.1 it 1.5 it 0.9 


1.0 - 


3.0 


(5.5 ± 1.1 ± 0.1) X 10-12 (2.2 ± 0.6 ± 0.2) x 10" 


-11 


1.6 it 0.3 it 0.03 


6.1 it 1.8 it 0.6 


3.0 - 


10. 


(6.2 ± 1.2 ± 0.1) X 10-" (2.2 ± 0.8 ± 0.2) X 10" 

Integral fluxitstatitsyst, cm— ^ s-^ 


-12 


1.9 it 0.4 it 0.03 


6.1 it 2.4 it 0.6 


>0.1 


(1.4 it 0.2t";4) X 10-^ (6.8 ± 0.7t°-'l) x 10~ 
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Note. — The S^Flux has been evaluated at the geometric mean of each bin. 



Table 3 

Integral flux (> 100 MeV) of the IC component 



Elongation 



Integral flux, 10 cm 



angle 


TS 


Elux±stat±syst 


Model 1 


Model 2 


Model 3 


0° - 5° 


93 


1.4 it 0.2^0 j 


1.15 


1.33 


1.78 


5° - 11° 


94 


2.5it0.3l°:i 


1.84 


1.97 


2.29 


11° - 20° 


43 


3.0it0.5lg;2 


3.11 


3.19 


3.50 


0° - 10° 


215 


3.7it0.4t°-j 


2.67 


2.95 


3.69 


0° - 20° 


259 


6.8 it 0.7l° j 


6.10 


6.48 


7.57 



Table 4 

Differential spectrum of the disk component 



Energ 


y range, 




Huxitstatitsyst, 




(i!;2Flux)±stat±syst, 


10^ 


MeV 


TS 


cm-2 MeV-i 




lO-'^ cm-2 MeV 


0.10 


-0.14 


59 


(3.1 it 0.5 it 0.6) X 10" 


-9 


4.3 ±0.7 ±0.9 


0.14 


-0.19 


68 


(2.4 ±0.3 it 0.5) X 10" 


-9 


6.1 ± 0.8 ± 1.3 


0.19 


- 0.27 


85 


(7.9 ± 1.1 it 1.6) X 10" 


10 


4.2 ±0.6 ±0.8 


0.27 


- 0.37 


100 


(3.2 it 0.4 it 0.6) X 10- 


10 


3.3 ±0.4 ±0.6 


0.37 


- 0.52 


222 


(2.2 it 0.2 it 0.4) X 10- 


10 


4.2 ±0.4 ±0.8 


0.52 


- 0.72 


288 


(1.2 it 0.1 it 0.2) X 10- 


10 


4.5 ±0.4 ±0.7 


0.72 


- 1.00 


243 


(5.2 it 0.6 it 1.0) X 10- 


11 


3.8 ±0.4 ±0.7 


1.00 


- 1.39 


282 


(2.8 it 0.3 it 0.6) X 10- 


11 


3.9 ±0.4 ±0.8 


1.39 


- 1.93 


295 


(1.6 ±0.2 ±0.3) X 10- 


11 


4.3 ±0.5 ±0.8 


1.93 


- 2.68 


193 


(7.5 ± 1.1 ± 1.5) X 10- 


12 


3.9 ±0.6 ±0.8 


2.68 


- 3.73 


167 


(3.4 ±0.6 ±0.7) X 10- 


12 


3.4 ±0.6 ±0.7 


3.73 


- 5.18 


143 


(2.0 ±0.4 ±0.4) X 10- 


12 


3.9 ±0.8 ±0.8 


5.18 


- 7.20 


85 


(8.2 ±2.1 ±1.6) X 10- 


13 


3.1 ±0.8 ±0.6 


7.20 


- 10.0 


21 


(2.2 ± 0.9 ± 0.4) X 10- 


13 


1.6 ±0.6 ±0.3 



Integral flux±stat±syst, cm ^ s 
>0.1 (4.6±0.2t,Y'j^) X 10-^ 



Note. — The B^Flux has been evaluated at the geometric mean of each bin. 
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Figure 3. Reconstructed spectrum of the background for the fake-Sun method (filled symbols, red) and for the simulated background sample (open symbols, 
blue) averaged over a 20° radius around the position of the Sun. 
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Residual counts 



Figure 4. Distribution of the residual counts: differences between tlie plioton counts in the map centered at the average fake-Sun and the high statistics count 
map derived from the model. A normal distribution with the same parameters is shown by the bold line. The map used to generate this distribution has a bin size 
of 0.3° centered on the solar position of the simulated sample. 
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Figure 5. Integral intensity profiles above 500 MeV for elongation angles <5° (top) and <20° (bottom). Points (red) are the observed counts, dash-dotted 
horizontal (black) line is the background, dotted (magenta) and dashed (green) lines are the point-like and extended components of the emission, con'espondingly. 
The solid (blue) line is the sum of the background and the two components of the emission. The shaded areas around the lines show total eiTor estimates. See 
text for details. 
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Figure 6. Electron spe ctrum used in the models of the IC emission compa red with the data from Fermi-LAT (small red filled circles, |Abdo et al.|2009bl 
[Ackermatm et al.|2010^ and AMS-01 (blue filled squares, [Alcaraz et al.|2000) . The thick solid (black) line is a fit to the Fermi-LAT electron spectrum at 1 AU ; 
this IS also the electron spectrum for r < 1 AU in Model i. the dash-dotted (black) and doted (black) lines show the electron spectra at at r = 0.3 AU in 
Models 1 and 2 calculated for <E>o = 400 MV, eq. pj and eq. |3|, correspondingly. The dashed (blue) line shows the demodulated local interstellar spectrum. The 
double-dot (magenta) line shows the local intersteflar spectrum of leptons as calculated by GALPROP i Ptuskin et al._2006j ; this spectrum is used in the model 
calculations of the diffuse Galactic emission. 
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Figure 7. Energy spectra of the IC emission for elongation angles <5° and <20° as observed by Fermi-LAT and compared with model predictions. Statistical 
error bars (larger) are shown in black; systematic errors (smaller) are red. 
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Figure 8. Intensity profile for the IC component vs. elongation angle compared with the model predictions. Statistical error bars (smaller) are shown in black; 
systematic en'ors (larger) are shown in red. To allow a direct comparison with the models, the model predictions are also shown binned with the same bin size as 
used for data. 




Figure 9. Energy spectrum for the disk emission a s obse rved by the Fenni-hAT. The curves show the range for the "nominal" (lower set, blue) and "naive" 
(upper set, green) model predictions by Seckel e t alT]^199T) for different assumptions about CR cascade development in the solar atmosphere (see text for details). 
The black dashed line is the power-law ttt to the data with index 2.11 it 0.73. 



